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I. INTRODUCTION 

As we approach its centenary [1], General Relativity 
(GR) remains our most promising theory for describing 
the gravitational interaction. It has had tremendous suc¬ 
cess in describing relativistic gravity in the solar system, 
and in binary pulsar systems [2]. When applied to cos¬ 
mology, it has resulted in the AGDM (or concordance) 
model of the Universe [3]. This model has the benefit 
of extreme simplicity, being based on a single universal 
Robertson-Walker geometry. It has also been extremely 
successful at reproducing a wide variety of cosmological 
observations, including supernovae [4], CMB anisotropies 
[5] , large-scale structure formation [6] , baryon oscillations 
[7], and weak leasing [8]. However, the ACDM model ap¬ 
pears to necessitate the largest fine-tuning problem that 
has ever occurred in physics, and the apparent coinci¬ 
dence that we live at the exact time when the cosmolog¬ 
ical constant is just starting to dominate the large-scale 
expansion [9]. These problems have led many authors to 
consider alternative theories of gravity [10], often based 
on modifications to the usual Einstein-Hilbert action that 
are expected to become large at low curvatures. 

One of the simplest ways to extend GR is to consider 
gravitational theories that have field equations that are 
derived from the metric variation of a Lagrangian den¬ 
sity that is a non-linear function of R: The so-called f{R) 
theories of gravity. Such theories have the pleasing prop¬ 
erty of being constructed entirely from the metric tensor, 
are well-motivated from attempts to improve the renor- 
malizability of gravity [11], and have often been consid¬ 
ered in early universe physics [12]. They have also been 
particularly well studied in recent years, as they poten¬ 
tially allow the late-time accelerating expansion of the 
universe to be interpreted as having a geometric origin, 
rather than being due to a violation of the energy condi¬ 
tions in the matter sector. 


Despite this wealth of work, very little effort has so 
far gone into understanding how weak-field systems can 
be embedded into cosmological solutions within this class 
of theories. Exact inhomogeneous cosmological solutions 
are very difficult to find, and Swiss-cheese-type construc¬ 
tions have only recently been considered [13]. In princi¬ 
ple, however, the cosmological behaviour should in some 
way determine the boundary conditions for the weak-held 
systems on small scales, and the cosmological expansion 
should be able to be viewed as being due to bodies on 
small scales falling away from each other. These ideas 
have been made explicit in General Relativity, where 
there exists a deep relationship between Newtonian force 
laws proportional to r~^ and r, and the dust-dominated 
and A-dominated Friedmann solutions [14]. It is the pur¬ 
pose of the current paper to consider the corresponding 
situation in f{R) theories of gravity that lead to acceler¬ 
ating expansion at late times. 

We begin by considering the Friedmann solutions of 
these theories, and demonstrate that there is great free¬ 
dom in the large-scale evolution of homogeneous and 
isotropic spaces, as long as one is allowed complete free¬ 
dom in choosing f{R). We then go on to consider how 
some particular theories, that have an evolution very 
close to AGDM today, generically exhibit growing os¬ 
cillations that lead to curvature singularities at finite 
times. In contrast, we then consider the cosmological 
models that emerge from joining small regions of per¬ 
turbed Minkowski space together using the appropriate 
junction conditions. The geometry inside each of the re¬ 
gions is expanded as a power series in the weak-field slow- 
motion limit. We find that the only theories that admit 
a viable weak-field limit within each of these regions are 
also forced to expand like ACDM on large scales. 

This result has several interesting implications. Firstly, 
it suggests that the freedom-in-principle that exists in 
the Friedmann solutions of f{R) theories of gravity does 
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not exist if one limits oneself to observationally viable 
theories. Secondly, it suggests that the potentially prob¬ 
lematic oscillations that can occur in their Friedmann so¬ 
lutions do not exist if the space-time can be described as 
slowly-varying weak-field perturbations about Minkowski 
space on small scales. Thirdly, it shows that viable f{R) 
theories of gravity cannot solve any of the problems as¬ 
sociated with the cosmological constant, or generalize its 
evolution in any way: Their large-scale evolution is iden¬ 
tical to ACDM, and the effective value of A must be 
constructed from parameters that appear in the gravi¬ 
tational action only. In this sense, there is only a nom¬ 
inal difference between late-time accelerating expansion 
in GR with A, and in viable f{R) theories. 

In this paper, we start in Section II by introducing 
the f{R) theories of gravity, as well as the mathematical 
tools that will be required build our models. In Section 
III we then continue to explore the freedom that exists 
in the FLRW models that are constructed within these 
theories, as well as the instabilities that can exist within 
an important class of them that is thought to be observa¬ 
tionally viable. In Section IV we then go on to consider 
cosmological models built from the bottom up. we find 
that the models that are built in this way, and that have 
an observationally viable weak-field limit, have severely 
constrained large-scale expansions. We summarize and 
conclude in Section V. 


II. PRELIMINARIES 


Let us now consider the tools required for our top-down 
and bottom-up constructions: That is, the Friedmann 
equations and the post-Newtonian expansion. 


A. Friedmann Equations 


The Friedmann equations are derived using the homo¬ 
geneous and isotropic Robertson-Walker geometry: 


ds^ = —dt^ + a^it) 


dr'^ 


1 — kr^ 


r'^{d9^+ sm'^ ed(j>^) 


( 3 ) 

where a{t) is the scale factor, and where k is the (con¬ 
stant) curvature of space. 

Substituting Eq. (3) into the field equations (2), and 
taking the matter source to be a perfect fluid, gives the 
Friedmann (constraint) equation as 


3H^=8ttp+pr, (4) 


where 

PR = - f) + ^ - 3Hf"R (5) 

is the conserved energy density of the effective f[R) dark 
energy, H = a/a is the Hubble parameter, and p and p 
are respectively the energy density and isotropic pressure 
of the perfect fluid (corresponding to standard matter 
fields). The Raychaudhuri (evolution) equation is 


The f{R) theories of gravity are a family of action- 
based theories, whose Lagrangian density is given by 

C=^fiR) + C^, ( 1 ) 

where R is the Ricci scalar, f{R) is a general differen¬ 
tiable (at least C^) function of the Ricci scalar, and Cm 
corresponds to the Lagrangian of standard matter fields. 
Here, and throughout, we use units in which c = G = 1. 

Variation of the action A = f d^x^J—gC with respect 
to the metric, yields the following gravitational field 
equations: 

- V,V^/' + W' = (2) 

where we used / = f{R), and where primes denote 
derivatives with respect to R. The energy-momentum 
tensor derived from varying the Lagrangian associated 
with standard matter fields is denoted by 

The Lagrangian (1) and field equations (2) are known 
to contain an additional scalar degree of freedom, that is 
not present in GR. This can be made explicit by confor¬ 
mal transformation of the metric [15, 16], or alternatively 
by a Legendre transformation of the theory. We will re¬ 
turn to this equivalency in Section V, but for the majority 
of the paper we will choose to keep the theory written as 
a function of the metric only, as in Eqs. (1) and (2). 


3H + 3H^ 


87r(p + 3p) + f - f'R (6) 

+3Hf"R -K 3f"R^ + 3 f'R . 


The energy density and pressure of the matter fields are 
related by the energy conservation equation: 

p + 3H{p + p) = 0, (7) 

while pr satisfies 

Pr + 3H{pr+pr) = 0 , ( 8 ) 


where 

pj, = f + 2Hf' + - f) . (9) 

The effective equation of state for f{R) dark energy is 
then given by 

— PR 

wr = — . (10) 

PR 

Here, and throughout what follows, we limit ourselves to 
the spatially flat geometries in which fc = 0. The other 
cases follow straightforwardly. 

The solutions to Eqs. (4)-(7) are the basis of the 
Friedmann-Lemaitre-Robertson-Walker (FLRW) cosmo¬ 
logical models, which will be discussed in Section HI. 
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B. Post-Newtonian Expansion 


The post-Newtonian expansion is a weak-field, slow- 
motion expansion about a Minkowski space background. 
It is the basis of almost all modern investigations of rela¬ 
tivistic theory in weakly-gravitating systems [2] , and has 
been shown to have an almost unreasonably wide field of 
applicability [17]. The expansion is based on assigning 
all quantities (geometric, kinetic and matter) orders of 
smallness in the parameter e = u/c, where v is the typical 
relative velocity of the matter fields that fill the space- 
time, and c is the speed of light. For the solar system, 
galaxies and galaxy clusters this parameter is expected 
to span the range 10“®-10“^. 

The geometry for the lowest order part of the post- 
Newtonian expansion can be written as 

ds^ = —(1 — 2(f>)dt^ -I- (1 -I- 2ip)dxidx^ + 2Aidtdx^, (11) 

where (/), ip and Ai are all gravitational potentials. The 
lowest order parts of the energy-momentum are then 
taken to be 


Too = P (12) 

Toi = -pvi (13) 

Tij = pStj , (14) 

where i and j run over spatial indices, and 6ij is the 
3-dimensional Kronecker delta. 

For ordinary matter content, the orders-of-smallness of 
the quantities in Eqs. (11)-(14) are given by 

, (15) 

and 

(16) 

(17) 

(18) 


A, 

P 

P 


Each and every time derivative is also taken to add an 
extra order-of-smallness to the quantity that it operates 
on, with respect to spatial derivatives, so that 


d/dt 

d/dx 


(19) 


These quantities are all dimensionless, making an assign¬ 
ment of an order of magnitude to them unambiguous. 
In the field equations, however, the terms that appear 
have dimensions of one-over-length-squared. The best 
one can do for such terms is therefore to say that they 
have an order-of-smallness when expressed in some choice 
of units for length. We choose these units such that spa¬ 
tial derivatives add no order-of-smallness to the objects 
on which they operate, which results in p ~ e^. This 
choice is for book-keeping purposes only, and does not 
affect any of the resulting equations. Any other choice of 
units would result in an identical set of equations. 


Of course, any additional fields that appear in the mod¬ 
ified theories of gravity, such an extra scalar in the case 
of f{R) theories, must also be expanded in powers of e, if 
the theory is to fit into the post-Newtonian framework. 
The procedure for solving Einstein’s equations (or the 
field equations of the theory in question, if different to 
GR) is then to solve each equation to lowest non-trivial 
order in e, before moving on to higher orders. 

The lowest orders in the diagonal components of Ein¬ 
stein’s equations are given to even orders in e, while the 
lowest order in the off-diagonal components are odd in 
e. The Newtonian limit of the equations of motion for 
massive particles is then given in terms of the O(e^) part 
of 4>, while the first post-Newtonian order is given by the 
square of the Newtonian terms, the O(e^) part of ip, the 
O(e^) part of Ai, and the 0(6"^) part of goo- The lowest 
non-trivial order in the equations of motion of massless 
test particles is given by the O(e^) parts of both (p and 
Ip, which are therefore sufficient to determine the post- 
Newtonian effects of gravity on rays of light. 

In Section IV we will construct a cosmological model 
by gluing together different regions of space-time, each of 
which is modelled using the post-Newtonian expansion 
described above. This will follow the methodology used 
in Ref. [14]. 


III. THE TOP-DOWN APPROACH: 
FLRW COSMOLOGIES 


What we call the ‘top-down approach’ to cosmology 
is, in fact, the usual way that cosmological models are 
constructed in the bulk of the literature. This approach 
starts by assuming that on sufficiently large scales the 
geometry of the Universe looks, in some sense, close to 
being spatially homogeneous and isotropic, so that it can 
be modelled using a Robertson-Walker geometry. It is 
then further assumed that this geometry should satisfy 
a given set of gravitational field equations, such as Eqs. 
(4)-(7). The former of these assumptions appears to be 
compatible with the majority of cosmological observa¬ 
tions, but the latter is a purely mathematical question. 
Proving its validity requires proving that a local theory 
of gravity, such as f{R) gravity, is applicable to a non¬ 
local averaged, or approximate, geometry. There are a 
number of extremely difficult technical questions that re¬ 
main to be answered in order to rigorously justify such 
an approach [18, 19]. Nevertheless, modelling the Uni¬ 
verse in this way appears to be able to account for the 
vast majority of cosmological observations with only two 
unknown matter fields. It is also an extremely simple 
way to model the Universe, and is used almost ubiqui¬ 
tously. lu this sectiou we will therefore briefly review the 
ELRW cosmological models of f{R) theories of gravity, 
emphasizing in particular the aspects that are of most 
relevance for the bottom-up constructions we will study 
in Section IV. 
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A. The Reconstruction Method 

This approach assumes that the expansion history of 
the universe is known exactly, and then inverts the field 
equations to deduce what classes of /(i?) theories give 
rise to this particular cosmological evolution [21]-[22]. It 
is useful for demonstrating the freedom that exists in the 
Friedmann solutions of f{R) theories. 

The reconstruction method works as follows. If we 
presume we know a = a{t), then we can take the usual 
expression for the Ricci curvature scalar, R = 6H+12H^, 
and invert it to find an expression for t as a function of 
R: 

t = gi{R). (20) 

This function can then be used to write all expressions, 
that were previously functions of t, as functions of R. 
The Friedmann equation (4) and conservation equation 
(7) then give 

3H[g^{R)]R[g,{R)]r + f " f) + { 

= Poa[5i(R)]-3(i+“), (21) 

where po is a constant, and where w = p/p has been 
taken to be the constant equation of state of the perfect 
fluid. This is a second-order differential equation that can 
be used to find the f[R) that produce the stated a{t). It 
works for almost any expansion history. Similar methods 
can be used to reconstruct f{R) from the functions a = 
( 72 (a), H = g^{H) and q = 94 ( 9 ) (where q = —aaja^ is 
the deceleration parameter) [ 22 ]. 

The freedom to produce such a wide array of differ¬ 
ent expansion histories has led to considerable interest in 
/(i?) gravity as a possible explanation of dark energy. In 
Section IV we will use bottom-up constructions to show 
that this freedom does not in fact exist, if the theory 
is required to simultaneously produce an observationally 
viable weak-field limit on small scales. All viable theories 
must expand like ACDM. 

B. Numerical Methods 

Numerical methods can also be used to investigate the 
Friedmann solutions of /(i?) theories. Of particular in¬ 
terest in this regard is the theory proposed by Hu & 
Sawicki [23]. This theory has proven to be extremely 
popular over the past few years, and is generated from 
the function 



where Rq^ b, d and n are all positive constants. The 
constant i?o has the dimensions of the Hubble parameter 


squared, and corresponds to the scale at which infrared 
modifications start to become important. Because both 
/' and f" are positive, the theory has a positive effec¬ 
tive gravitational constant and does not initially appear 
to posses any obvious instabilities. On the surface, this 
theory (and others like it [24, 25]) therefore appear to 
represent viable descriptions of a cosmic acceleration at 
late-times, driven by nothing more than the curvature of 
space-time. 

An example of the numerical solutions to Eqs. (4)-(7) 
for the Hu & Sawicki theory (22) with n = 3 is presented 
in Fig. 1. In producing these plots we have assumed that 
H{z) and q{z) take on their ACDM values at the present 
day {z = 0), and have chosen Rob/d = GDAiJQ. While the 
Hubble parameter as a function of redshift is very close 
to the ACDM value, it is clear that growing oscillations 
are present in the deceleration parameter, and that these 
oscillations translate into an effective equation of state 
for dark energy that varies wildly. In fact, it turns out 
that these oscillations lead to curvature singularities at 
finite times. This instability was pointed out by Frolov 
[26], and has since been studied numerically [27, 28]. 

These oscillations represent a severe drawback for the 
viability of such theories. Their suppression requires a 
careful fine-tuning of the parameters of the model, on 
top of what is required to set the value of the effective 
cosmological constant. On the other hand, they poten¬ 
tially provide a way of testing these /(i?) theories against 
supernovae and growth rate measurements [29]. A fur¬ 
ther difficulty for the n = I case, with the same initial 
conditions as above, is that the matter dominated epoch 
at high redshifts does not even exist [30]. 

In Section IV we will show that these oscillations do 
not occur at leading order in the bottom-up construc¬ 
tions, as long as the slow-motion weak-field limit holds on 
small scales. This shows that, for perfectly homogeneous 
and isotropic FLRW models, the post-Newtonian limit 
must be destroyed on small scales when the oscillations 
start. However, it also suggests the intriguing possibility 
that the growth of structure in the real Universe could 
suppress the unstable oscillations. This is because non¬ 
linear structures keep the Ricci curvature scalar (which 
is a single-valued locally-defined quantity) in the high- 
curvature regime at all points in space. A corollary is 
then that the expansion histories of universes that obey 
such theories must be constrained to be even closer to the 
ACDM model than had previously been thought. We will 
return to these points in Section V. 


IV. BOTTOM-UP CONSTRUCTIONS: 

A PATCHWORK UNIVERSE 

What we mean by ‘bottom-up constructions’ are cos¬ 
mological models in which the geometry of space-time in 
the vicinity of astrophysical objects is the thing that is 
solved for first, and given primary importance. This can 
usually be done up to the application of suitable bound- 
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Figure 1: The Hubble, deceleration and effective equation of 

state of dark energy wn defined by equation (10), for the Hu &; 
Sawicki model with n = 3. Here b = 1, d = 0.1 and we have taken 
= 0.76. The condition Rgb/d = gives Rq = 0.456i7Q. 


ary conditions, for the partial differential equations in¬ 
volved. Such boundary conditions are imposed by patch¬ 
ing together a large number of sub-horizon sized regions 
of space-time using the Israel junction conditions. These 
same junction conditions simultaneously specify the mo¬ 
tion of the boundaries between regions, and hence can 
be used to extrapolate the large-scale cosmological ex¬ 
pansion of the space-time as a whole. In this approach 
to cosmological modelling, the field equations need only 
be applied locally (i.e., not to any non-locally averaged 
geometries), and the cosmological expansion can there¬ 
fore be viewed as an emergent phenomenon, as it should 
be. The approach we will use here was originally pro¬ 
posed for GR in [14], and was further developed in [31]. 
It was adapted to the study of some f{R) theories of 
gravity in [13]. Here we will further adapt it to study the 
emergence of accelerating expansion in f{R) gravity. 

The situations we will primarily be concerned with 
are ones in which matter is clumped into highly non¬ 
linear structures on small scales. We will then take 
a sub-horizon sized region of space, which we will call 
a “cell”. We have in mind that these cells should be 
polyhedral. We will treat the geometry of space-time 
inside each of these cells as being well approximated 
by the post-Newtonian expansion, as outlined in Sec¬ 
tion IIB. We will then consider many such cells, all the 
same shape, that are invariant under discrete rotations, 
and that are patched together about reflection symmet¬ 
ric boundaries. On very large scales the global geome¬ 
try of space-time can then be thought of as statistically 
homogeneous and isotropic, while being very inhomoge¬ 
neous and anisotropic on small scales. At no point in this 
process do we assume a global Robertson-Walker back¬ 
ground. 

While this approach is quite flexible, for definiteness we 
could consider each cell to contain a galaxy-sized amount 
of mass (~ 10^^ solar masses), and to have a spatial 
extent that is similar to the inter-galactic separation in 
the actual Universe (~ 1 Mpc). If we then assume that 
the rate of expansion of our model Universe is of a similar 
order of magnitude to that of the real Universe (~ 70 km 
s“^ Mpc“^), then at the edge of each cell we will have a 
Hubble flow given by 

u ~ 2 X 10"^c. (23) 

The typical gravitational potential for structures of this 
size, separated by this distance is then 

^ ~ 5 X 10-® . (24) 

We therefore have 



where e ~ 10“"^. This is sufficient to be able to apply 
a weak-field and slow-motion post-Newtonian expansion 
of the field equations. A similar process could be consid¬ 
ered for clusters of galaxies, which as long as the average 
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density is unchanged, would result in an expansion pa¬ 
rameter of order e ^ 10“^. 

Crucially, when considering general f{R) theories of 
gravity, we will not be able to make any assumptions 
about the smallness (or not), of / and /'. We will also 
have to refrain from expanding these quantities around 
i? = 0, even though our perturbative expansion within 
each cell will be around Minkowski space. This is done 
so as to include as many f{R) theories as possible, some 
of which contain parameters that are small, and some of 
which are only expected to admit convergent power series 
expansions around finite values of R. We will, however, 
have to assume that / and f can be expanded in our 
order-of-smallness parameter e (if this is not true then 
the theories cannot be treated perturbatively at all). In 
this case we have 

f = fio) ^ f{2) ^ (26) 

and 

f = fio) + f'i2) + 0{e^), (27) 

where superscripts in brackets denote the order of small¬ 
ness of a term in the e expansion. As usual in post- 
Newtonian expansions, we have expanded these scalars 
in even powers of the expansion parameter e. The con¬ 
vergence of these expansions can be verified a posteriori. 

The junction conditions we will use when joining cells 
together will be those found by Deruelle, Sasaki and 
Sendouda [32] in the absence of surface layers, and when 
/VO and /" ^ 0: 


[labt 

= 0 

(28) 

[Kabt 

= 0 

(29) 

[Rt 

= 0 

(30) 

{dyf]t 

= 0, 

(31) 


where "fab and Kab are the first and second fundamental 
forms on the boundary, where a,b... are indices on the 
1 -|- 2-dimensional boundaries, and where dy is a normal 
derivative on the boundary. 

For the lattice construction described above these junc¬ 
tion conditions provide Neumann boundary conditions 
for the field equations within the interior of each cell. 
This can be contrasted with Swiss cheese constructions, 
in which the boundary conditions imposed on these 
equations would be of Cauchy-type, and the imposition 
of “cosmological boundary conditions”, which are often 
taken to be of Dirichlet-type (with the idea that weak- 
field solutions slowly turn into cosmological solutions at 
very large distances). 

The geometry that will be used will be given by Eq. 
(11). Initially, it may seem that the zero-order equations 
here are completely trivial (i.e., that they are immedi¬ 
ately solved as long as Minkowski space is a solution of 
the theory). However, we should bear in mind that many 
of the theories of interest to us at present are those that 


have been constructed to have interesting behaviour in 
f when R is small. We therefore cannot immediately 
assume anything about the smallness of the lowest order 
non-zero parts of / and /', but instead must determine 
their magnitude through the field equations (2). 

From the zero-order part of the field equations it can 
be seen from + ri? = 0, = 0 and + 

Tzz^ = 0 that we must have 

= ffJ = 0. (32) 

The equations Txy'^ = = 0 also give 

/•'(O) _ fl{0) — f'(0) _ Q ('QQ'l 

J ,xy J ,xz ,yz 

Substituting this into the Tq® = 0 equation then gives 

= 0 , (34) 

so that the leading order term in / must be O(e^) or 
smaller. The leading order term in /', on the other hand, 
must be given by the integral of Eqs. (32) and (33): 

/'(o) = a{t) + (3i{t)x + I32{t)y + , (35) 

where /3i, (32, and a are arbitrary functions of t. How¬ 
ever, it can immediately be seen from the junction con¬ 
dition [dyf]^ = 0, and the imposition of reflection sym¬ 
metric boundaries, that we must have dyf' = 0 at every 
point on the edge of every cell. Using Eq. (35), this 
means that we must have l3i = P 2 = Ps = 0. 

The leading order behaviour of f must therefore be 
spatially homogeneous, so that 

f = a{t) + 0{e^). (36) 

This result suggests that we can consider two separate 
classes of theory: 

• Class I: Theories that can be written as /(i?) = 
aR + F{R), such that /' = a + F', where a is a 
constant (not a function of R). 

• Class II: Theories with R = Rq + 0{e^), where 

Ro = Ro{t) ~ 

In Class I we have all theories in which is not a 
function of R, but is instead constructed from constant 
parameters that appear in the gravitational Lagrangian. 
In this case we can integrate to find f{R) = aR + F{R), 
such that F F' . The Class I solutions contain 
most of the viable f{R) theories that have recently been 
considered in the literature (e.g. Hu & Sawicki [23], Ap¬ 
pleby & Battye [24], and Starobinsky [25]). 

In Class H we have all theories in which f'Fi is a non¬ 
trivial function of the Ricci curvature scalar, R. If this is 
true, together with the fact that /' is a function of t only, 
then the lowest order part of R must also be a function of 
t only. In this case we can write R = Rq + 0(e"^), where 
Rq = i?o(i) ^ Any space dependent contributions 
to R^‘^'> would not be compatible with the supposition 
a{t) = a{R). 

Let us now consider each of these two classes individ¬ 
ually. 
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A. Class I Theories 

In this case the lowest-order parts of F and F' are both 
O(e^). The equations Tqq^ = Sirp and = —Sirp then 
give, to lowest order, 

aV^ip + ^V^F' = -Anp+^F (37) 

aV^cj) - = -47rp - ^F, (38) 

where denotes a spatial Laplacian operator. Likewise, 
the Txy^ = TxV = Ty’^z = 0 equations give 



= 0 

(39) 

^,xz ^'4^,xz ^^,xz 

= 0 

(40) 

^,yz ^'4^,yz 0^0,yz 

= 0 . 

(41) 

These latter equations require 



F' = «((),- 7 /)) -p ^ (ci(t,x) + C 2 (f, 2 /) +C 3 (t, 2 ;)), 

(42) 


where ci, C 2 and C 3 are (so far) arbitrary functions of t, 
with X, y and z, respectively. Combining Eqs. (37), (38) 
and (42) gives 

F = V2(ci+C2 + C3), (43) 

which in turn allows us to note that, because F^x = 
F'R^x ^ we also have at lowest order that 

Cl,XXX — 0; C2,yyy = 0 aud ^‘i^ZZZ ~ 0- (^4) 

This tells us that ci, C 2 and C 3 all contain terms up 
to quadratic order in their respective spatial arguments 
only. Equation (43) then becomes 

F = Fo(t), (45) 

and allows Eqs. (37) and (38) to be written as 

“ —47rp+ -Fq (46) 

V2 (^acj> - ^F'^ = -Airp - ^Fq. (47) 

At this point it is convenient to consider the = 
—Sirpvi equations: 

-aV'^Ai - -aAjji + atp^i + -F'j = Airyvi. (48) 

Substituting this into the time-derivative of Eq. (46) 
then gives 

47r(puj),i = -471/9 + ^Fo, (49) 

which on using the Eulerian continuity equation of hy¬ 
drodynamics for a perfect, non-viscous fluid gives 

Fq = constant. 


Equations (46) and (47) can now be seen to have the 
solutions 


i’ 


1 

a 

1 

a 



(C + Ir'). (51) 


where U and U denote the usual Newtonian potentials, 
including cosmological constant, and are given implicitly 
by the following expressions: 
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V^U = — 47 rp-|--Fo 

(52) 

2 - 1 

V^U = -4^p--Fo. 

(53) 


Now, the large-scale expansion of the cosmological 
model obtained from patching together our weak-field re¬ 
gions is given by [14]. The outward acceleration of each 
point on the boundary of every cell is 

X = dy4'+0{e‘^). (54) 


Or, defining the average outward expansion by 

f xdA 


(55) 


where the integrals here are over the surface of a cell face, 
we can write 


d _ 47rGeff M ^ 
a 3 720^ 3 


(56) 


where Ges = 372/nQ;7i, where M = J pdV is the integral 
of the energy density over the volume of the cell, where 
n is the number of faces to a cell, and where 

1 _ GeffFg /c:y\ 

=-^— (57) 

is the contribution to the cosmological constant from the 
F{R) part of the action. We have also defined the area 
of a cell face and the volume of the cell to be A = 710 ^ 
and V = 720 ^, where 71 and 72 are constants. Of course, 
if the cell maintains its shape during the evolution, as 
expected, then we have d = x. Notice that both Geff and 
A/j are constants, to the required order. 

The parts of Eq. (56) that look like dust and a cos¬ 
mological constant come from the normal derivative of 
U on the boundary of the cell. The terms involving the 
normal derivative of F', however, do not contribute to 
the large-scale evolution at lowest order at all, because 
the junction condition [dyf']~^ = 0 evaluated at O(e^) 
gives dyF' = 0. The bottom-up cosmological solutions of 
all theories in this class therefore behave like GR plus A, 
with corrections appearing at O(e^) only. 

The only new behaviour possible in this sub-class arises 
in the interior of each cell, and is contributed by the 





possible presence of an O(e^) term in F'. Any such term 
must obey the trace of the field equations: 

= ^R+^Fo-^p. (58) 

The solutions to this equation may or may not obey a 
screening mechanism, but either way do not affect the 
leading-order behaviour of the large-scale expansion. 

As Fq behaves as the effective cosmological constant, 
and because Fq must be constructed from constants that 
appear in the gravitational action only, this class of the¬ 
ories does not appear to solve any of the problems asso¬ 
ciated with the cosmological constant (i.e., why should 
Fq cancel with the vacuum energy that appears on the 
right-hand side of field equations, and why should the 
value of Fq be of the same order of magnitude as the 
energy density in matter at the present time?). Beyond 
this, one may also notice that there are no signs of the 
oscillations that appear in the FLRW solutions to these 
theories. 


B. Class II Theories 


In this class of theories we are unable to directly inte¬ 
grate /' in order to find the lowest-order part of f{R). 
We therefore resort to writing 

/ = /(2)+o(e4) (59) 

and 

/' = a + /'(2)+0(e4)^ (60) 

where a = a(t) ~ 1. In this case, the equations Tqq^ = 
Sirp and 2Tqq^ -|- = Sirp give 

1/(2) _ v2//(2) _ ^^2^ ^ 

and 

-- 3d 4-- 4aV2V’= Sttp , (62) 


(21 

and the T-b ' = 0 equations can be integrated to find 
/'(2) =a{(j)-'ip) + ^{ci{t,x) + C2it,y)+C3{t,z)), (63) 


where Ci, C2 and Cq are arbitrary functions of their ar¬ 
guments. These equations can then be manipulated to 
obtain 


aS/^ip = —■ 


3 ^ e-i 3 

+ ^V^(ci +C2 + C3) 


and 


= - 




a 


-^V^(ci -bC2 -bca). 


(64) 


From the definition of R we also have 

aR = 2a{V^(j3-2V‘^ip) 

= /(2) + 2d-V2(ci+C2 + C3). (66) 

Now, because we have in this class of theories that R^"^^ = 
R^'^\t), we must also have that (as this is 

the leading order term in /, which is a function of R 
only). Equation (66) then shows that ci, C2 and C3 must 
be at most quadratic in x, y and z, respectively, so that 
V^(ci -b C2 -b C 3 ) is a function of t only. 

Using this information, Eqs. (64) and (65) can be writ¬ 
ten 


Stt 

aV^'0 = ——p + homogeneous terms (67) 

O 

and 

aV^(j) = --^-p + homogeneous terms , (68) 

O 

which have the solutions 
2U 

ijj = - -b cosmological terms (69) 

30: 

and 


. 41/ , . , 

(p = -b cosmological terms 

3a 


(70) 


where U is the Newtonian potential that satisfies 

= -47rp, 


and where “cosmological terms” refers to potentials that 
are independent of p, and proportional to r^. It can 
immediately be seen that all theories in this class have 
the post-Newtonian parameter 

1 

7PPN = ■^ = 2 ’ 


where cosmological terms have been neglected. The 
Chameleon Mechanism cannot apply to this class of the¬ 
ories, as it requires R to be inhomogeneous for it to be 
effective. This can be seen from the trace of the field 
equations, which can be written as 


V2//(2) ^ (72) 

Stt 

= —Td + homogeneous terms . 

O 


The solution for is then simply 


/d2) = 


2U 


cosmological terms, 


(73) 


( 65 ) 


which does not permit the type of non-trivial distri¬ 
butions that are required to implement the Chameleon 
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mechanism, and suppress the propagation of the scalar 
degree of freedom. 

The fact that 7 ppn = 1/2 is also evident directly from 
the definition of R, which gives 

R = 2VV-4VV (74) 

327r 

= —-—(1 — 27 ppn)p + homogeneous terms . 

3a 

For this equation to be simultaneously consistent with 
R = R{t), and an inhomogeneous p = p{t,x,y, z), we 
clearly require 7ppn = 1/2. 

We therefore find that theories in this class, while po¬ 
tentially admitting new and interesting Friedmann so¬ 
lutions, cannot be considered observationally viable, as 
they conflict with experimental probes of post-Newtonian 
gravity, which currently require [33] 


7PPN - 1 = (-1.7±4.5) X 10"^ (75) 

This is inconsistent with yppN = 1/2 at more than lOOOcr. 

V. DISCUSSION 

We have studied the ability of f{R) theories of grav¬ 
ity to account for the present day accelerated expansion 
of the Universe using the top-down and bottom-up ap¬ 
proaches to cosmological modelling. The top-down ap¬ 
proach admits a wide variety of possible evolutions. They 
also tend to admit oscillations about background expan¬ 
sions that are close to that of the ACDM model. These 
oscillations are a potential way to distinguish between 
f{R) cosmology and the ACDM model, but they also 
appear to signify an instability within the FLRW models 
of f{R) theories. 

Within the bottom-up constructions we identified two 
possible classes of theory, which we named Class I and 
Class II. The Class I theories are potentially observation- 
ally viable, and include most of the well-studied theories 
that have recently appeared in the literature [23]- [25]. 
The expansion history that emerges within this class of 
theories is indistinguishable from the ACDM model, and 
is not plagued by the potentially problematic oscillations 
that are known to occur in the Friedmann solutions of 
these theories. This appears to be a type of strong back- 
reaction, and signifies that either the growth of structure 
suppresses the oscillations that occur in the homogeneous 
and isotropic solutions, or that the occurrence of oscilla¬ 
tions destroys the e expansion that we have used to de¬ 
scribe the gravitational field on small scales. This could 
mean that either the fluctuations become large or rapidly 
changing, so that the e expansion is no longer convergent, 
or that the Minkowski space background becomes unsuit¬ 
able. In either case the weak-field slow-motion treatment 
of gravitational fields must be said to be no longer ap¬ 
plicable. If this is the case then it is extremely difficult 
to see how the theory could be considered to be observa¬ 
tionally viable. 


The Class II theories have the potential to describe new 
and interesting large-scale behaviour, due to the new cos¬ 
mological terms in Eqs. (69) and (70). However, it can be 
seen that they must all have a weak-field limit that is not 
compatible with observations of post-Newtonian gravity, 
as they have 7 ppn = 1/2. The Chameleon mechanism 
cannot apply to this class of theories, as the lowest-order 
part of R must be a function of time only. So, while the¬ 
ories in this class are potentially capable of displaying a 
variety of cosmological expansion histories, we are led to 
conclude that they are not observationally viable on all 
scales. 

Our key result is therefore the following: All f{R) the¬ 
ories that are potentially compatible with observations 
have expansion histories in the late Universe that are in¬ 
distinguishable from GR with A. Moreover, the value of 
the effective cosmological constant in these models must 
be a combination of the bare cosmological constant, the 
vacuum energy density (pvac), and a possible combination 
of other constant parameters that appear in the gravita¬ 
tional action of the f{R) theory in question: 

47r 

AefFective — A^are Pvac “b • (76) 

a 

The fine-tuning that is required to get AeSective at the 
level required in order to explain observations is there¬ 
fore not at all alleviated by generalising the gravitational 
Lagrangian from the Einstein-Hilbert one to the f{R) 
family. To say that theories in this class cause accelerat¬ 
ing expansion would also appear to be misleading: The 
only parts of the gravitational Lagrangian that are con¬ 
tributing to the acceleration at late times are those parts 
that, to leading order, are constants. 

A natural question that arises is whether the results 
found above also apply to scalar-tensor theories of grav¬ 
ity, as the f{R) theories are known to be equivalent to 
a class of these theories after a Legendre transformation. 
To be more explicit, if we make the definition 

^ = f, (77) 

and if we assume ip = p{R) is invertible, so that we can 
make the further definition 

= \ [77(7^)75 - /(7?(</j))] , (78) 

then we are led to the following gravitational Lagrangian: 

C=^[pR-2V{p)]. (79) 

This is the sub-class of scalar-tensor theories that have 
vanishing coupling constant, a; = 0. In this case the field 
equations can be written 


pG^y -I- {V^Vap + V)g^u - y^NuP = , (80) 

and the scalar field propagation equation becomes 

(81) 




2 dV 
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Taking the trace of Eq. (80), and using Eq. (81), then 
gives 


R = 




(82) 


This is an algebraic relationship between the Ricci cur¬ 
vature scalar of the space-time, and the scalar field in the 
gravitational theory. It shows that in scalar-tensor the¬ 
ories with to = 0 and V(ip) ^ constant that if we know 
information about the curvature of the space-time (by 
assuming it can be described by a weak-field expansion 
on small scales, for example), then we also know some¬ 
thing about the scalar degree of freedom in the theory. 


This was a key point used in finding the results above 
for /(i?) theories of gravity, and so we also expect our 
results to apply to the a; = 0 sub-class of scalar-tensor 
theories. This does not, however, mean that they should 
be expected to apply to scalar-tensor theories in general, 
as there is in general no algebraic relation between R and 

tf. 
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